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SUMMARY 

This  investigation  was  performed  at  Princeton  University  in  the 
period  from  February,  1951  to  Jime,  1951* 

The  original  intention  was  an  experimental  investigation  of  the 
feasibility  of  xislng  a  powered  model  in  harnessed  circling  flight  to  evalu- 
ate the  longitudinal  stability  of  a  f\ill-scale  airplane.  The  airplane 
selected  vas  the  Cessna  1^0  tvo-place  light  airplane,  since  its  stability 
derivatives  and  other  parameters  had  been  determined  and  vere  at  hand.  A 
one-tenth  scale  flying  model  vas  constructed. 

The  method  of  attacK  vas  to  have  been  as  follows: 

A.  Determination  of  stick-fixed  neutral  point,  N^  frcm  steady- 
state  flight  tests* 

B.  Szperlmental  determination  of  the  model's  frequency  response- 
T  by  obtaining  a  transient  response  of  9  to  an  elevator  pulse  iiiput 

and  converting  this  to  the  frequency  response  by  the  method  outlined  in 
Reference  11. 

C.  Correlation  of  A  and  B  vith  theory. 

No  quantitative  data  vere  obtained  for  A  and  B  due  to  mechanical 
difficulties  and  lack  of  time.  S\ifficient  qualitative  data  vere  obtained, 
however,  to  show  agreement  betveen  the  model's  transient  response  aixd  that 
predicted  by  theory;  to  point  up  several  difficulties  peculiar  to  this  type 
of  testing;  and  to  indicate  recosDamendatioiis  for  further  development  of  the 
method. 

It  vas  foiind  that  the  scaJ^  chosen  gave  too  small  a  model  from  co]>- 
siderations  of  required  mechanical  Installations  and  dynamic  similarity. 
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Transient  response  tests  of  relatively  short  duration  were  found 
to  be  the  only  feasible  type  due  to  the  presence  of  vake  disturbances. 
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IHTRODUCTIOH 

With  the  increasing  nioiber  of  nev  and  reidically  different  configura- 
tions of  airplanes  being  introduced  for  sonic,  transonic  and  convertaplane 
purposes  it  seems  that  seme  method  of  evaluating  the  lov-speed  longitudinal 
stability  characteristics  in  unaccelerated  forward  flight  is  desirable  prior 
to  fixing  a  design.  This  could  most  easily  be  done  in  model  form.  A  tech- 
nique such  as  this  vould  allow  the  inertia  forces  on  the  model  to  be  inter- 
preted to  full  scale  results  and  vould  also  be  valuable  in  determining  gust 
loadings. 

The  Air  Materiel  Command  at  Wright  Patterson  AFB  established  a  test 
facility  which  gave  promise  of  yielding  such  data. 

An  analytical  study  (Beference  3)  was  made  by  AMC.  From  this  it 
was  concluded  that  certain  of  the  longitudinal  stability  characteristics 
could  be  evaluated,  and  that  even  where  the  dynamic  stability  of  the  p>vt;goid 
mode  was  greatly  modified  by  the  harnessing  forces  on  the  model,  the  model 
response  could  be  related  to  fxill  scale  response  by  certain  mathematical 
operations.  Unfortunately  the  test  facility  was  decomnissioned  before  any 
data  were  taken  to  substantiate  the  above  conclusions,  and  at  the  time  of 
ccomencement  of  this  investigation,  had  not  been  re-established. 

This  facility  made  use  of  powered  models  flown  in  harnessed  circling 
flight;  a  type  of  flight  known  as  "control-line  flying"  which  is  vised  as  a 
sport  or  hobby  by  certain  of  the  general  public. 

It  was  the  pxirpose,  therefore,  of  the  present  investigation  to 
provide  a  limited  euaount  of  information  on  which  to  evaluate  harnessed  cir- 
cling flight  of  models  as  a  testing  method,  and  to  establish,  if  possible, 
testing  techniques  which  will  complement  the  method.  To  this  end  a  one-tenth 
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scale  model  of  the  Cessna  1^0  light  two-place  airplane  constructed.  Reasons 
dictating  choice  of  this  airplane  were  that  complete  static  flight  tests  and 
vlnd  tunnel  tests  of  the  Cessna  140  have  been  previously  made  here  by  the 
Aeronautical  Engineering  Department  of  Princeton  University  and  therefore 
complete  data  vere  available.  This  investigation  vas  limited  to  ccoiparing 
the  neutral  pointy  as  obtained  from  level  flight  testing  of  the  model,  with 
the  neutral  point  from  flight  test  of  the  fxill-scale  airplane;  and  comparing 
the  frequency  response  of  the  model's  pitch  angle  to  that  obtained  from  the 
equations  of  motion  of  the  airplane .  The  experimental  frequency  response 
was  to  be  obtained  from  the  analysis  of  the  transient  response  to  a  triangu- 
lar pulse  of  elevator  deflection. 

These  limitations  were  dictated  by  lack  of  facilities,  funds  and 
time. 

The  data  were  to  be  obtained  from  static  stability  flight  tests  of 
the  model  with  three  center-of -gravity  positions  and  from  transient  response 
tests  at  one  center-of -gravity  position  {21, df)  mac). 

It  Is  believed  that  this  investigation  is  the  only  one  of  its  kind 
to  date.  Although  the  final  results  may  not  be  conclusive  due  to  rather 
primitive  equipment,  it  may  indicate  a  method  and  provide  a  i>art  of  the  in- 
formation necessary  for  a  complete  analysis  of  the  method,  which  of  course, 
requires  numerous  tests  on  sftany  more  of  the  longitudinal  stability  aspects. 
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STKBOLS  AHD  ASSUMPTIONS 

AB     -    Wiog  aspect  ratio 
AC     -    Aerodynamic  center  of  the  wing 
^     _    MCu'l  -  Slope  of  the  lift  cvtyq  of 

the  wing  (two  dinenslonal) 


- = m 


(per  degree) 

-  Slope  of  the  lift  curve  of  the 
w 

wing  (three  dimensional) 

(per  degree) 


a*  ■=-(^^\  -  Slope  of  the  lift  curve  of 


the  tail  (three  di&ensional) 
(per  degree) 
D  -  Wing  span  —  feet 

C  -  Mean  Aerodynaiaic  Chord  of  the  Wing  —  feet 
C.G.-  Center  of  Gravity 

Ow  ~  u  oyxc  -  Wire  drag  coefficient 


Cp  =   "^  ^  P w  -  Drag  Coefficient 

Cp-  Effective  drag  coefficient  {Cq-Cj) 

^T  ~  v^  '  Thrast  coefficient 

C^  -  Parasite  drag  coefficient 


D 


-f 


Cjj-  Induced  drag  coefficient 
C^    -  ^^  -  Lift  coefficient 

^^t  -  1=1  -  Lift  coefficient  of  the  tail 

C,v^  z   ^  -  Pitching  Moment  coefficient 

of  the  airplane 
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Cj,   r  ir^  -  Change  in  drag  coefficient 

with  change  in  angle  of 
attack. 
Cd*  =   *^^o  -  Change  in  effective  drag  coefficient 

vith  change  in  angle  of  attack. 
/^  (f  -  Angle  between  flying  vires  and  datim  plane 
^  ^«  ~  f-£.'-  "  Change  in  lift  coefficient 


vith  change  in  angle  of 
attack  —  per  radian 
*^«  ~  - — *"  -  Change  in  pitching  moment 

coefficient  vith  change  in 
eu3gle  of  attack  —  per  radian 
^m  1  =  ^^>^  -  Change  in  pitching  moment 

VXt/  )  coefficient  vith  rate  of  change 

of  angle  of  attack  vith  respect 
to  ^/y 

•^de  ~  — •         -  Change  in  pitching  moment  co- 

^^    r)  efficient  vith  rate  of  change 

of  pitch  angle  vith  respect 
to  */^ 
r         ^     dC  "  ^^1""^®  i^  pitching  moment  co- 

^  i>  efficient  vith  change  in 

elevator  angle. 
D   -  Total  drag  of  the  airplane  —  Ihs 
d   -  Differential  operator 
e  -  0.869  -  Airplane  efficiency  factor 


PRINCETON  UNIVERSITy  PAGE  7 

AERONAUTICAL  ENGINEERING  LABORATORY  REPORT    i8i 


F       -  Centrifxjgal  force  in  level  flight 

r  F  'T 

r  -  -r  -— -r,  -  Hon-dimens  ional  line  force 

R  pvS 

2   -  Acceleration  di*e  to  gravity  —  32,2  ft/secS 

~      -  -  Hon-dimena ional  inertia  parameter 

-  Airplane  moment  of  inertia  about 
2         -5  the  Y  axis  —  slug  ft2 

1  or  j     -    \^r 

mac   -  Mean  aerodynamic  chord  —  per  cent 

k^  -  Radius  of  gyration  about  the  Y  axis  —  ft 

L  -  Total  lift  of  the  airplane  —  lbs 

Ji^   -  Tail  length  (airplane  CG  to  AC  horizontal  tail)   —  feet 

r\    -  Moment  about  the  Y  axis  —  lb -ft 

y^    -  Mass  of  airplane  —  slugs 

"^^P*^  -  Miles  per  hour 

Nq  -  Stick-fixed  neutral  point 

n   -  Normal  acceleration  (along  the  7.  axis  at  the  CG)  —  "G"  units 

P  -  Period  of  oscillation  —  seconds 

^  -  /^  p  V  ^  -  Dynamic  pressure 

R   -  Radius  of  flight  circle. 

S   -  Wing  area  —  ft^ 

Sju  -  Horizontal  tail  area  —  ft^ 

Tv  -  Time  to  damp  to  l/2  amplitude  —  sec 

X    -  Time  —  sec 

-^  -  Speed  ratio  ~1    ,  incremental  change  in  airspeed  divided 

V 
by  airspeed 

V  -  True  airspeed  —  ft /sec 
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W  =  5*JL*  -  Volumetric  coefficient  of  tlxe  tail 

S  c 

W  -  Gross  weight  (average)  —  pounds 

gRSEK  SYMBOLS 

^  -  Angle  of  attack  of  ving  —  degrees 

o<^  -  Angle  of  attack  of  the  horizontal  tail  —  degrees 

^^   -  Elevator  deflection  —  degrees 

€.  -  Angle  of  dovzxirash  —  degrees 
'l)t  -  ^  -  Airplane  tail  coefficient 

6  -  Angle  of  pitch  (angle  between  X  axis  aM  horizontal  —  degrees) 

A  -  Boot  of  stability  quartic 

/vv\ 
/^  -  — 5 —  -  Relative  airplane  density 

factor  (non-dimensional) 

P   -  Mass  density  of  air  —  slugs/ft^ 

T  -  — r-rr  -  Time  parameter  --  seconds 

(^  "     Phase  angle  »  degrees 

CO  -  Angular  frequency  --  radians/sec 

CO  -  Angular  frequency,  radians  per  X/ 

A  variable  with  svibscript  zero  indicates  the  maximum  value  of  that 
variable  during  axi  oscillation  and  may  be  either  positive  or  negative. 

Vertical  bars  indicate  the  conventional  absolute  value  of  any  variable. 

SIGH  COHVEimOH 
Reference  Axes 

Normal  stability  axes  (vind  axes)  are  used. 
A  -  axis  -  longitvidinal  axis  -  positive  forward 
Y-  axis  -  lateral  axis  -  positive  along  right  wing 
Z-axis  -  vertical  axis  -  positive  down 
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Llne^  Displacements 

A  llaear  displacement  stlong  a  positive  reference  axis  is  considered 
positive, 

Angxilar  Displacements 

An  angular  displacement  which  is  clockvise  vhen  viewed  frcm  the 
origin  looking  along  a  positive  reference  axis  is  considered  positive. 

Velocities  and  Accelerations 

Velocities  and  acceleration,  either  linear  or  angular,  are  con- 
sidered positive  in  the  same  sense  eus  the  corresponding  displacements. 

Control  Surface  Deflections 

Positive  elevator  angle  is  associated  with  a  dovmrard  movement  of 
the  elevator  trailing  edge. 

Phase  Angles 

A  negative  phase  angle  signifies  a  time  lag  of  the  response  to 
the  forcing  function. 

A  positive  phase  angle  signifies  a  time  lead  of  the  response  to 
the  forcing  function. 

ASSnMPTIOWS 

In  the  theoretical  analysis  and  in  the  evaluation  of  the  stability 
derivatives  the  folloviog  assumptions  are  made: 

1.  Stability  aixes  (vind  axes)  are  xised. 

w 

2.  The  airplane  is  initially  in  level  flight.   (This  xobbjob   for 
stability  axes   cos6=i,  stt^B-Q    ) 

3>  All  disturbances  are  small,  so  that  the  equations  of  motion  may 
be  considered  linear. 


L 
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k.     The  lift  caused  by  the  horizontal  tail  is  negligible. 
^.  Stoall  lateral  disturbances  vhich  may  exist  do  not  affect 

the  aerodynamic  forces  or  pitching  mcaaents  in  the  plane  of 

symmetry  • 
6.  Forces  due  to  pitching  velocities  are  negligible. 
7*  Cheuige  in  pitching  mcanent  due  to  the  rate  of  change  of 

elerator  angle  is  negligible. 

8.  Only  longitiidinal  motion  is  considered.  The  moment  of  inertia 
is  that  about  the  Y  axis. 

9.  Criteria  for  dynamic  similarity  were  based  on  equal  C  ^    of 
model  and  airplane  (Beferences  ^  and  10). 

10.  Standard  model  flight  conditions  used  in  the  solution  of  the 
equations  of  xaotion  vere: 

V  =  58.7  1^ 
B  =   37.5  I* 
at  datum  level  =  J*  ft.  above  the  floor. 
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DSSCRIPTIOH  OF  THE  MODEL  AMD  INSTRIMENTATIQH 

The  ttLodel  used  In  this  inveetigatlon  vas  a  ozie-tenth  scale  model 
of  tlie  Cessna  1^  tvo-place  light  airplane.  It  vas  constructed  by  the 
authors.  The  constrtiction  consisted  primarily  of  balsa-wood,  i)aper-covered 
and  doped;  a  construction  similar  to  spin  tunnel  test  models  (Beference  9), 
Original  intentions  had  been  to  meike  the  model  dynnmically  similar  on  the 
basis  of  equal  lift  coefficient  (Beferences  3  and  10).  Hovever,  the  choice 
of  airplane  and  the  scale  used  soon  izidicated  that  this  vo\ald  be  physically 
iBq;>ossible  with  available  equipoent,  due  to  the  lov  veight  requirement  and 
small  size  of  the  model.  The  choice  of  airplane  vas  fixed  by  considerations 
previously  mentioned,  and  scale  vas  fized  by  (1)  budgetaxy  reasons  plus 
(2)  the  requirement  of  an  indoor  flying  site  to  reduce  stray  vind  disturbances 
on  the  model.  It  vas  decided  to  stay  as  close  as  possible  to  the  veight 
limitation  and  measure  the  resulting  moment  of  ineirbia  of  the  model,  using 
that  and  the  actual  veight  in  the  non-dimensiooal  equations  of  motion  of  the 
airplane  of  the  harnessed  model.  This  modification  in  no  vay  affects  the 
ability  to  Judge  the  effectiveness  of  this  type  of  flight  as  a  test  method. 

The  model  vas  povered  by  a  Bantam  .19,  cubic  inch  displacement, 
l/lO  HP  (@  9000  EPM),  internal  cooxbustion,  air-cooled,  model  airplane  engine 
vith  glov  ignition  (Figure  I  ).  This  ignition  saved  veight  by  requiring  no 
ignition  components  in  the  model;  merely  a  l^^-volt  battery  as  a  starting 
device.  The  engine  vas  movuated  inverted  to  preserve  the  contours  of  the 
upper  covling  and  it  projected  out  of  the  lover  cowling  a  short  distance 
(Figure  2  ).  A  pressiare  type  fuel  feed  syst«a  of  ccamnercial  make  vas  used 
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to  prevent  variation  of  engine  speed  due  to  change  of  fuel  head  in  flight 
(Figiire  I  ). 

In  order  to  give  acese  to  the  mechanism,  engiz^  and  friel  system 
the  vlng,  as  veil  as  varioxis  hatches  vere  made  removable  (Figure  ^  ) . 

In  order  to  clarify  the  operation  of  the  elevator  pulsing  mechan- 
ism it  is  necessary  to  describe  the  control  of  the  model  in  flight.  In 
normal  flight  the  model  flies  in  a  circle  around  the  operator  on  two  vire 

lines  which  lead  froaa  a  bellcrank  (Figure  ^  )  in  the  fuselage,  through  a  guide 

To  ft 
near  the  wing  tip  handle  in  the  hand  of  the  operator.  The  model  is  flown 

with  some  yaw  to  the  outside  of  the  circle  which  is  produced  by  offsetting 

the  rudder  and  by  attaching  the  bellcrank  behind  the  C.G.  The  euaount  of 

this  yaw  is  also  controlled  by  adjxisting  the  line  guide  on  the  wing.  A  push 

<l^vaVor  horn 

rod  from  the  bellcrank  links  with  tKe  so  that  a  motion  of  the  operator  tend- 
ing to  increase  tension  in  the  top  line  from  the  handle  and  decrease  tension 
in  the  bottoaai  line  is  transmitted  as  an  up  deflection  of  the  elevator. 

The  requirements  on  the  elevator  for  the  pvilse  input  are  that  it 
(1)  deflects  from  the  trim  elevator  position,  (2)  returns  to  this  position, 
and  (3)  remains  locked  during  the  ensuing  transient  motion.  To  this  end  a 
third  wire  from  the  operator  to  the  model  was  employed.  When  a  level  equili- 
brium flight  condition  was  attained,  the  operator  transferred  tension  from 
the  two  flying  wires  to  the  third  wire.  The  point  of  attachment  of  this 
wire  was  slightly  aft  of  the  model's  C.G.  in  order  to  maintain  slight  outward 
yaw  and  hence  line  tension.  The  first  result  of  tension  on  the  third  line 
was  the  application  of  a  braJse  on  the  elevator  push  rod  (Figure  M   ).  This 
locked  the  elevator  in  the  equilibrium  position  previously  established.  As 
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the  tension  vas  increased  to  coimteract  the  full  centrifijgal  force  of  the 
models  the  third  line  also  triggered  a  z*ubber-povered  pulsing  mechanism 
(Figure  ^  ).  The  cam  in  this  mechanism  displaced  the  aft  portion  of  the 
push  rod,  working  against  a  spring  vhich  returned  the  elevator  to  the  original 
locked  position.  The  cam  vas  c\xt  to  give  a  pulse  of  trianguleir  shape  vith 
time,  having  a  duration  of  approximately  .0^  seconds.  By  adjusting  the  x>o8i- 
tion  of  the  cam  fo Hover  on  the  push  rod,  the  amplitude  of  the  pulse  vas 
varied.  The  cam  contacting  the  follover  closed  an  electrical  circuit  throu^ 
a  third  vire  vhich  fired  a  fla^h  bulb.  After  the  elevator  had  been  pulsed, 
the  mechanism  vas  autcanatically  re-cocked  for  another  run  vhen  the  tension 
on  the  third  line  vas  released  and  the  elevator  a^in  controlled  by  the  oper- 
ator. 

Measvirement  of  elevator  angle  vas  necessaiy  in  both  phases  of  the 
testing.  In  the  steady-state  tests  the  angle  required  for  eqvillibriimi  at 
each  velocity  had  to  be  known,  and  in  the  dynamic  tests  the  determination  of 
the  exact  size  of  the  pulse  vas  necessary.  This  angle  vas  to  be  measured 
vith  a  Giannini  Micro-torque.  This  instrtuoent  is  an  electrical  potentiometer 
eaa^loying  a  rotating  contact  of  very  low  mechanical  torque  gradient.  Angular 
positions  are  measxired  electrically  by  suitable  calibration.  In  this  applica- 
tion the  microtorque,  acting  as  a  variable  resistance,  vas  mounted  in  the 
fuselage  and  attached  by  a  linkage  to  the  elevator  (7ig\are  H   ).  Recording 
elevator  position  required  tvp  leads  from  the  model  to  the  center  of  the  cir- 
cle. Enameled  steel  flying  vires  vere  used  for  this  purpose.  From  the  oper- 
ator's position  the  leads  vere  continued  to  the  recording  apparatxas  on  the 
ground  (Figure  b  ).  The  steady-state  elevator  position  vas  recorded  by 
calibrating  a  voltmeter  to  elevator  position.  A  recording  oscillograph  vas 
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employed  to  record  the  pulse. 

Ho  In-flight  meosxirements  vere  possible  in  practice  for  tvo  reasons. 
First,  the  vibration  induced  by  the  znmning  engine  vas  sufficient  to  destroy 
the  electrical  continiiity  of  the  micro-torq;ja .  Second,  the  friction  intro- 
duced into  the  elevator  system  by  the  zaicro-torque  linkage  made  control  of 
the  Biodel  marginal. 

Weights  adjustable  fore  and  aft  vere  used  to  alter  the  model  C.Q. 
position  in  running  the  steady-state  phase.  Hormal  C.6.  for  the  model  cor- 
responded to  noznal  C.G.  for  the  full  scale  Cessna  and  vas  at  27. d^  mac. 
This  position  vas  used  for  the  dynamic  phases. 

Dynamic  pitch-angle  data  vas  recorded  by  means  of  a  motion  picture 
camera  tracking  the  model  from  the  center  of  the  circle.  This  camera  vas 
mounted  on  a  heavy  trix>od  eq\iix>ped  vith  leveling  bolts  and  a  revolving  pljaten. 
A  flash  bulb  in  the  3^6.  line  circviit,  placed  so  that  its  flash  vould  shov  on 
the  film,  vas  used  to  mark  the  instant  of  application  of  the  pulse.  The 
camera  vas  equipped  vith  a  cross  hair  reticule.  With  the  camera  field  verti- 
cal and  the  camera  adjusted  to  track  in  the  horizontal  azimuth  plane,  these 
cross  hairs  vere  a  reference  for  both  pitch  angle  measurement  and  vertical 
displacement.  "Die  camera  speed  vas  calibrated  for  a  time  reference  during 
the  transient  response. 

The  use  of  aji  armory  vas  obtained  for  a  testing  site.  This  reduced 
the  amount  of  stray  air  currents  as  compared  vith  an  outdoor  site,  and  Improved 
the  chances  of  obtaining  reproducible  results. 
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TBST  PROCEDURE 

The  CO.  location  of  the  model  was  found  by  suapending  it  from 
three  positions  having  vide  angiilar  separation  about  the  pitch  axis.  The 
C.G.  vas  at  the  intersection  of  plumb  lines  hung  from  the  x>oint8  of  sus- 
pension. 

The  model's  moment  of  inertia  about  the  pitch  axis  was  determined 
by  swinging  it  as  a  coopound  pendulvsa.  The  method  is  described  in  Reference 
6  and  calculations  given  in  Appendix  G. 

This  C.G.  location  i21»&f>  mac)  and  moment  of  inertia  (.012  slug  ft^) 
were  used  in  calculating  parameters  for  the  equations  of  motion  of  the  model. 

Ihie  to  restrictions  imposed  by  conditions  at  the  flying  site  (low 
ceiling  lamps)  and  by  the  light  weight  of  the  models  it  was  necessary  to  fly 
it  froaa  a  37  •5  foot  radius  instead  of  the  50  foot  radi\is  originally  contem- 
plated. T^ils  resulted  in  a  higher  rate  of  yaw,  but  enabled  the  outward  yaw 
angle  of  the  model  as  set  by  the  line  guide  on  the  wing^  to  be  reduced  be- 
cause of  the  increase  In  centrifugal  force  with  decrease  in  radius.  The 
shorter  radius  was  objectionable  from  the  viewpoint  of  slipstream  disturbeoices 
as  the  wake  had  less  time  to  die  out  before  being  flown  through  again  by  the 
model. 

A.  Static  Pheise 

No  static-phase  runs  were  made  due  to  the  above-mentioned  mechanical 
difficulties.  For  the  determination  of  the  stick-fixed  neutral  point  the 
model  was  to  have  been  flown  in  steady  level  flight  at  a  given  C.G.  position 
and  velocity.  The  velocity  would  have  been  measured  by  timing  a  niaaber  of 
laps.  The  trim  elevator  position  would  have  been  measured  with  the  calibrated 
voltmeter.  A  number  of  flights  were  to  be  made  at  this  C.G.  position  with 
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various  speeds.  The  speed  would  have  been  varied  by  a  combination  of  engine 
adJustBient  and  propeller  change.  This  procedure  of  flying  at  various  speeds 
vas  to  be  repeated  for  a  number  of  C.G.  locations. 

B.  Dynamic  Phase 

With  the  pulsing  mechanism  wound  vip,  the  model  was  flown  level  emd 
the  speed  was  measured.  The  electrical  leads  were  connected  to  the  opera- 
tor's handle  from  the  recording  apparatus,  and  the  recording  oscillograph 
stairted.  The  camera  was  euLso  started,  emd  tracked  the  model  under  manual 
control.  The  pulse  was  triggered,  the  pitch  angle  data  being  recorded  by 
the  camera  and  the  elevator  pxalse  by  the  oscillograph.  The  record  was  con- 
tinued until  the  stable  transient  died  out. 

KBSULTS 

No  static  phase  data  were  possible  because  the  micro-torque  angle 
indicator  would  not  record  under  the  vibrating  conditions  produced  by  the 
engine  of  the  model. 

A  total  of  fifteen  flights  were  made  during  which  the  model  was 
adjusted  for  controllable  flight  and  qualitative  data  were  observed  for  the 
dynamic  phase.  These  consisted  of  approximate  measurement  of  speed  and 
visual  observation  of  the  following:  characteristics  of  the  phugoid  mode, 
response  to  5,  10,  and  25-degree  elevator  pulses,  and  the  behavior  of  the 
wake  made  visible  by  exhaust  smoke. 

All  elevator-locked  flights  were  above  datum.  Speeds  observed 
were  between  30-50  mph.  The  characteristics  of  the  phugoid  were  erratic. 
Its  period  was  observed  to  vary  between  1.5*2.5  seconds.  In  some  instances 
it  was  seen  to  diverge  sharply.  In  one  instance,  during  kO   seconds  of 
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elevator-locked  flight,  several  phtigold  oscillations  vere  seen  to  commence 
and  damp  out,  under  the  Influence  of  vake  dleturbancea  or  other  effects. 

The  2^-degree  elevator  pulse  vas  found  necessary  for  a  visually 
observahle  pulse. 

At  the  last  pulse  during  one  of  the  runs,  the  energy  remaining  In 
the  mechanism  vas  lov  enough  to  appreciably  Increase  the  duration  of  the  pulse. 
The  result  vas  a  greatly  Increased  transient  leading  to  a  violent  phugold 
oscillation  vhlch  eOmost  terminated  the  entire  experimental  Investigation. 

DISCUSSION 
Static  Phage: 

Various  power  settings  are  reqiilred  during  tests  for  elevator  pover 
and  stick-fixed  neutral  point  using  harnessed  flight.  This  results  in  changes 
in  elevator  pover  and  the  necessity  for  evaluating  thrust  outputs.  Thus  data 
reduction  is  contplez,  and  since  these  characteristics  may  be  found  more  easily 
from  vind  tunnel  tests  it  vas  concluded  that  their  calcxxLation  frcm  harnessed 
flight  tests  is  inadvisable.  Another  deterrant  is  external  disturbance  due  to 
flight  in  the  slipstream  of  the  previous  lap. 

Dynamic  Phase; 

The  initial  intentions  of  this  phase  vere  to  collect  aiid  analyze 

the  flight  data  using  the  photographs  of  the  time  history  of  the  pitch  angle 

to  plot   r-   against  time.  This  plot  vas  then  to  be  graphically  analyzed  to 
o 

find  the  frequency  response  as  described  in  Reference  11.  This  experimentally 
determined  frequency  response  curve  vas  to  be  ccmpared  vith  the  frequency 
response  ttcex  the  analysis  of  the  transfer  function  of  the  equations  of  motion 
for  the  conditions  of  harnessed  flight  and  of  free  flight.  Since  the  equations 
of  motion  developed  in  Appendix  A  are  non-dimensional  and  Include  the  harnessing 
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forces  as  a  separate  term,  they  apply  equally  veil  to  a  model  or  to  an  air- 
plane dynamically  similar.  Because  the  small  size  of  the  model  prevented 
the  attainment  of  a  dynamically  similar  mcment  of  Inertia,  the  equations,  as 
solved,  refer  not  to  the  full  scale  Cessna  l^i-O  but  to  an  airplane  having  all 
the  Cessna  longitudinal  stability  derivatives  and  characteristics  (Eeferences 
7  and  8)  except  for  the  mass  and  pitching  moment  of  inertia.  These  two  char- 
acteristics, found  by  direct  measxirement  of  the  model  and  calculated,  are 
those  of  an  airplane  dynamically  similar  to  the  model  bviilt  for  this  investi- 
gation. 

It  was  assimied  from  the  Information  in  the  AMC  report  (Reference  5) 
that  the  stabilizing  and  daioplng  effect  of  the  harnessing  line  forces  on  the 
phugoid    oscillation  would  eliminate  this  mode  from  any  consideration  in 
the  analysis  of  flight  data  and  woxild  prevent  it  from  becoming  a  problem  in 
the  testing  technique.  This  did  not  prove  to  be  the  case.  Although  on 
flights  with  the  elevators  under  control  of  the  operator  the  model  appeared 
stable  with  the  CO.  at  27«93t  mac;  when  the  elevators  were  locked  the  model 
exhibited  a  very  short  period  phugoid  (1^  to  2?  seconds  period)  when  flown 
at  elevations  of  about  10'  above  datum,  Fiirthermore  this  phugoid  oscilla- 
tion was  qviite  unstable,  increasing  to  double  amplitude  in  less  than  one  per^ 
lod  on  niaierous  occasions.   It  was  sufficiently  unstable  to  make  pxilsing 
the  elevator  at  the  datum  level  (4*  above  the  floor)  when  flying  at  kO  miles 
per  hour  possible  only  at  the  risk  of  desti*oying  the  model.   It  was  decided 
that  to  enable  data  to  be  taken  safely  a  higher  datum  plane  was  necessary, 
i.e.,  the  operator  standing  on  an  elevated  platform  while  flying  the  model. 
See  sketch  on  page  19.  Since  the  location  permitting  this  was  not  available, 
no  quantitative  data  were  obtained  at  present  writing.  Instead,  the  amphasls 
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of  the  inY«Btlgation  wae  shifted  to  detenmination  of  the  caiiaes  of  the  erron- 
eoua  prediction  auid  to  analyzing  methoda  which  vould  make  the  fulfillment  of 
the  alms  of  the  original  Investigation  possible. 


Po>.^o>.irv\    plArve 


5ketch  show/inQ  metKod  3f  Ya\5iKva 
dcxtawv  plana 


In  an  att^qpt  to  find  the  causes  for  the  discrepancy  between  the 
observed  flights  and  the  theory,  the  assxaoptions  used  in  writing  the  equa- 
tions of  motion  (Appendix  A)  were  re-examined.  It  was  confirmed  that  the 
angle  of  the  lines  leaving  the  model  in  the  drag  direction  could  be  neglected. 
This  angle  was  less  than  2^  \jnder  the  most  extreaae   flight  conditions.  It 
was  also  found  that  the  angle  in  the  gravity  direction  due  to  the  weight  of 
the  lines  was  also  negligible.  Hence  the  assumption  that  the  lines  extend 
in  a  straight  line  frcm  the  operator  to  the  model  is  valid  even  in  so  light 
a  model.  This  is  diw  primarily  to  the  short  radius  employed.  The  conditions 
that  the  model  must  be  charged  with  3 A  ^^  ^^^   drag  of  the  hfiLmessing  wires 
and  that  the  mass  of  the  wires  and  their  inertia  effects  upon  non-steady 
motions  can  be  neglected  (due  to  the  small  mass  of  the  wire)  were  also  veri- 
fied. Calculations  supporting  the  preceding  are  shown  in  Appendix  D. 
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Eovever,  observatlonji  shaved,  that  at  the   short  radius  flovn  (which  was 
necessitated  tj   indoor  floor  space  and  celling  obstructions)  and  at  the  scale 
speed,  the  assxai^tlon  that  the  circling  flight  would  produce  no  aerodynamic 
effects  on  the  longlt\idlnal  motion  was  not  completely  Justified.  Calcula- 
tions show  that  the  speed  ratio  between  the  velocity  at  the  outer  wing  tip 
and  the  inner  wing  tip  is  1.1?^  for  the  model  used  in  this  investigation. 
For  cotBoparlson,  to  show  that  the  short  radius  used  here  was  no  more  detrimental 
to  this  condition  than  the  radius  used  at  the  Wright  Field  testing  facility, 
a  calculation  was  made  based  upon  their  B-29  model  (Beference  ^).  This  re- 
sulted in  a  speed  ratio  of  1.22  between  outer  and  inner  wing  tip.  Therefore 
the  model  used  in  the  present  investigation  was  less  affected  by  this  varia- 
tion of  velocity  with  spaxrwlse  position  than  the  large  one  flown  at  Wright 
Field  on  I50'  lines.  This  speed  ratio  effect  would  produce  a  very  cosiqBlicate d 
solution  to  the  equations  of  motion.  Therefore  in  the  solution  of  the  equations 
of  motion  it  is  omitted,  with  the  reservation,  however,  that  at  least  portions 
of  the  di8crex>ancles  between  the  theoretical  solution  for  the  frequency  response 
and  the  experimental  must  be  chargeable  to  this  amission.  An  additional,  aero- 
dynamic effect  produced  by  the  circling  flight  path  is  the  Interference  effect 
of  the  model's  own  wake.  Flight  observations  showed  an  intermittent  disturb- 
ance which  was  most  easily  noticed  by  a  sudden  inward  banking  oscillation  of 
the  model.  This  was  determined  to  occur  when  the  model's  outer  wing  passed 
through  the  slipstream  remaining  from  a  previous  circuit.  The  banking  oscilla- 
tion was  produced  by  the  fact  that  in  counterclockwise  circling  flight  the  out- 
board wing  was  being  influenced  by  an  iipward  component  of  both  the  slipstream 
and  the  inboard  wing  tip  vortex.  Confirmation  of  this  reasoning  was  assisted 
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by  the  fact  that  the  exhaust  emoke  of  the  engine  made  the  sllpetream  vake 
visible.  This  vake  vas  seen  to  leave  the  model  approximately  teoigent  to 
the  flight  circle  vlth  some  dovxnrash.  See  sketch  belov. 


Ration  0  + 


t^o<^«!  a1  Acj  on 


Skttch  s,hi>wiv\^  wc»^ke  effect^ 


Since,  at  scale  speed  the  model  passes  the  same  azimuth  point  on  the  flight 
circle,  and  the  vake  appeared  to  viden  as  fast  or  faster  than  the  curvature 
of  the  flight  circle,  the  model  crosses  this  turbulent  area  unless  it  is 
flovn  above  the  previous  path  by  a  reasonable  distance. 

All  of  the  foregoing  has  not  indicated  any  conclxaalve  reason  for 
the  vmstable  phugoid  vhich  vas  observed.  It  vas  noted,  hovever,  that  in 
both  the  Einalytical  solution  made  by  Braun  in  Reference  5  and  in  the  initial 
solution  for  the  present  investigation,  the  familiar  a8SunQ>tion  of  gliding 
flight  condition,  namely  C»p  «  0  ,  vas  used.  The  solution  obtained  for  this 
condition  (Appendix  C)  using  the  equations  developed  for  harnessed  flight 
developed  in  Appendix  A,  gave  as  roots  for  the  modes  of  motion: 
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A.=  -1.438,   A^=  -5.35 
^3.M=  -  •001  i  .2631 
This  indicates  a  very  lightly  damped,  stable  phiigoid,  vith  a  period  of  ^.1 
seconds  vhlch  dampe   to  I/2  amplitude  in  30  cycles.  This,  of  course,  was 
coangpletely  contraury  to  experimental  observations.  It  is  8J.80  of  interest 
to  note  that  the  short  period  mode  has  degenerated  into  two  real  roots. 
This  indicates  that  the  model  was  overdamped,  caused  by  too  low  a  moment  of 
inertia.  This  moment  of  inertia  discrepancy  had  been  expected.  (See  Des- 
cription of  Model). 

In  order  to  evaluate  the  effect  of  thrust  on  the  equations,  the 
characteristic  equations,  both  with  and  without  the  harnessing  force,  were 
re-solved  using  the  following  thrust  estimate:  Cj  »  Cjj  for  all  flight  con- 
ditions close  enough  to  the  equilibrium  condition  such  that  the  small  angle 
approximations  hold.  This  assizraption  makes  Ct)»s  0  and  C*  =  0.  Solution  of 
the  characteristic  eq\xation  (Appendix  C)  yielded  the  following  roots: 

A,::  -1.705         A^=  -5.88l^ 
A^^^  r  .0055  ±  .265  jL 

This  solution  indicates  a  trend  in  the  right  direction  since  it  gives  an 
unstable  phugoid  of  essentially  similar  period  to  the  glide  condition  (5.06 
seconds),  but  which  increases  to  double  its  amplitude  in  5*32  cycles.  This 
corroborates  the  observed  instability  of  the  phugoid  although  it  does  not 
show  the  rapid  divergence  which  was  observed  in  tests.  It  also  shows  that 
the  power  has  only  a  very  slight  effect  on  the  short  period  mode. 

To  check  the  validity  of  using  the  previous  thrxist  assunrption  an 
additional  solution  was  made  of  the  characteristic  eqviation  (Appendix  C) 
assuming  that  C^  was  constant  at  the  value  of  C  at  scale  speed.  This  makes 
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C-  «  0  and  C-.       -  Cp  .  The  changes  produced  in  the  coefficients  of  the 
characteristic  equation  between  this  and  the  previous  assumption  vere  so 
!nn»^''T  that  they  vere  within  the  accuracy  of  the  calcxilation  method  xtsed  to 
factor  the  characteristic  equation.  It  was  therefore  concluded  that  the 
previous  assvonption  was  sxifficiently  precise  for  the  size  of  the  factors  ajid 
for  the  method  used  in  extracting  the  roots. 

As  a  further  correlation  of  the  period  of  the  phugoid  obtained 
from  the  equations  with  the  period  from  observations  the  value  of  the  period 
was  corrected  for  flight  at  conditions  other  than  at  datum  level.  This  was 
accomplished  by  the  approximate  relationship   co  ^  .  /  bii-  zz     ^t  from 

V  a    s/2r 

Reference  1.  Using  the  standard  flight  conditions  and  simple  trigonometry,  . 
a  force  of  -.119  l^s  per  foot  of  height  from  datum  was  calculated,  (measuring 
directions  and  forces  in  airplane  stability  axes  direction).  This  is  con- 
verted into  a  C  -  correction  factor  by  dividing  by  the  weight  and  adding 
it 

one:  C,^  •  Cj^  (l  -V-  .O78  ^   )•  Therefore  denoting  by  subscript  (  )2^  condi- 
tions other  than  at  datim,  the  following  approximate  relationship  can  be 
expressed: 

^'  -  ^  -  C   -     Cu      -      * 

P      ^'    Cl,    Cu(»4-.078M      1  4-  .078W 

From  the  above  it  can  be  seen  that  flying  the  model  at  13  fe»t  above  datum 
approximately  halves  the  period.  Therefore  the  observed  results  agree 
reasonably  with  the  analysis  insofar  as  the  period  of  the  phugoid  is  concerned. 
Also  the  quality  of  the  stability  of  this  mode  is  in  agreement,  namely  they 
are  both  \instable.  The  only  discrepancy  is  the  rapidly  divergent  oscillation 
which  was  observed  but  which  failed  to  materialize  from  the  theory.  This  is 
unexplained. 
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The  foregoing  dlscuasion  has  been  primarily  about  the  phugoid 
mode  Bince  this  mode  indicated  the  discrex>ancy  in  the  initial  assimrptions 
upon  which  the  equations  of  motion  were  solved.  Nov  that  reasonable  agree- 
ment has  been  established  between  the  theory  and  the  qualitative  obseirva- 
tions,  the  a88UDQ)tions  which  produced  this  agreement,  namely  C^  ^  0,  Cq^'  ^  t 
will  be  used  in  sill  further  analysis.  The  phugoid  mode  is  in  itself  of  no 
interest  in  this  investigation  except  as  a  coo^licating  factor,  since  it  is 
most  markedly  affected  by  the  harnessing  force  of  the  lines  (Appendix  C). 
Examination  of  the  short  period  roots,  on  the  other  hand,  shov  practically 
no  difference  between  those  of  harnessed  and  unharnessed  flight.  It  is  this 
fact  that  appears  to  make  this  testing  method  desirable  for  evaluation  of 
dynamic  stability  above  the  frequency  of  the  phugoid* 

It  is  cozisidered  insufficient  to  compare  the  dynamics  of  ccoplex 

systems  on  the  basis  of  the  characteristic  equations  alone.  Such  an  analysis 

oixly  considers  the  natural  modes  of  motion  of  the  system  when  disturbed,  but 

not  under  the  action  of  a  forcing  function.  Since  the  normal  dynamic  testing 

methods  involve  the  application  of  a  forcii^  function,  it  is  more  reasonable 

to  ccopare  harnessed  and  free  flight  on  the  basis  of  the  frequency  response. 

This  firequency  response  can  be  obtained  in  the  conventional  manner  by  solving 

the  transfoxmed  equations  of  motion  for  the  transfer  function  —   as  shown 

o 

in  Appendix  B.  The  transfer  functions  to  a  step  input  were  reduced  to  the 
freqxienoy  domain  by  the  graphical  logarithmic  method  of  Reference  2,  and 
Illustrated  in  Appendix  K.  The   freqxiency  response,  as  obtained  from  the 
analysis  of  the  equations  of  motion  for  both  harnessed  and  unharnessed  flight 
are  plotted  on  Figures  7  ^md    S     .  These  are  the  curves  with  which  the 
experimental  frequency  response  was  to  have  been  compared.. 
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Examination  of  the  frequency  response  curves  Figures  7  and  6 
show  that  above  a  freqiiency  of  1.0  radian  per  ^/^  the  frequency  response 
for  harnessed  flight  is  identical  in  shape  with  that  of  the  free  flight,  but 
the  amplitude  of  the  former  is  less.  The  fact  that  this  difference  in  anqpli- 
tude  is  a  constant  factor  is  quite  evident  frcan  the  logarithmic  plot  of 
Fig\ire  10  .  On  this  plot  there  is  a  constant  2  decibel  gain  difference. 
Converting  this  gain  difference  into  amplitude  it  was  determined  that  the 
amplitude  of  the  harnessed  model  vas  80?t  of  the  correct  free  flight  valxie  at 
all  frequencies  above  co  =  l.o.  This  loss  of  amplitude  is  not  unreasonable 
since  the  lines  offer  restraining  forces  which  would  tend  to  limit  the  ampli- 
tude of  oscillation  in  response  to  a  forcing  function. 

Since  the  aforementioned  loss  in  eaaplitude  is  tied  up  in  numerous 
terms  of  the  transfer  function,  it  would  be  difficult  to  evaluate  it  directly 
frcan  the  harnessing  force.  However,  it  would  appear  that  the  change  in  am- 
plitude is  a  fimction  only  of  the  mass,  speed,  emd  flying  radius  of  the  model. 
IRxerefore  it  should  be  possible  to  evaluate  this  percentage  change  for  a  var- 
iety of  permutations  of  these  items  and  tabulate  the  result  to  be  used  in  ad- 
justing eacperimentally  determined  frequency  response  curves  to  free  flight. 

In  order  to  see  if  the  transfer  function  of  the  harnessed  flight 
Z     gave  a  transient  response  which  was  similar,  qualitatively,  to  that 

o 

observed  when  the  model  elevator  was  pulsed;  and  also  to  establish  desirable 
frame  speed  for  the  recording  camera,  it  was  considered  necessary  to  solve 
for  the  transient  response.  In  order  to  siuqplify  this  caJ.culation,  a  square 
piilse  was  assumed  insteeid  of  the  triangular  one  actually  imposed  on  the 
model.  It  shoiild  be  noted  that  a  triangular  pxilse  was  chosen  for  the  model 
to  enable  the  simplified  method  of  Reference  11  to  be  used  in  converting 
from  the  exi>erimental  transient  response  back  to  the  frequency  response. 
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The  mechanical  details  of  detemdnlDg  the  transient  response  from 

the  transfer  f\mctlon  exe  covered  In  tlie  literature.     Numerous  methods  may 

be  used.     The  one  chosen  for  use  here  Is  the  Seavlslde  expansion  vhlch  Is 

described  In  References  1  and  4,  namely, 

Given  the  transfer  function       -(A^     —    LiA) 

'^       &(M 

and  roots  of  G(X)  are  K,     X*>X3*''''^n  noDB  of  which  are  zero  or  re- 
peated roots, 

--  1^*^  --  '£  ^Ixf^.^  ^'-  ^'^^^^  =  '^^ 

This  solution  is  the  time  response  to  a  unit  step  input  (Appendix 
E).  By  using  the  method  of  Reference  k   the  response  to  a  p\ilse  of  finite 
vidth  can  be  found.  This  method  considers  the  pulse  as  being  made  i:^  of  a 
positive  unit  step  followed  at  time  At  by  a  negative  xmit  step.  By  s\:g>er- 
posing  these  two  inputs  a  unit  pulse  of  width  At  is  obtained.  By  the  prin- 
ciple of  the  Convolution  Theorem  the  output  is  then  the  result  of  si:g>erposiiB 
the  outputs  to  these  step  inputs.  These  outputs  are  identical  in  shape  but 
of  opposite  sign  and  separated  in  tima  by  an  interval  Ait.  .  These  operations 
were  accoanplished  on  logarithmic  coordinates  and  the  results  are  plotted  on 
Figures  i  I  emA   I  2.  .   These  are  time  histories  of  the  pitch  angle  9  per 
unit  elevator  input.  The  solutions  are  for  a  pulse  of  .Oh   seconds,  .10 
seconds,  and  .IH   seconds.  They  show  that  the  pitch  changes  of  the  tremsient 
are  of  reasonable  amplitude  for  these  two  pulse  widths,  provided  that  leurge 
elevator  deflections  are  used.  Ihls  corroborates  the  vis\3ally  observed  ef- 
fects. They  eulso  show  that  the  pulse  must  be  kept  short  or  the  transient 
response  due  to  the  short  period  mode  is  negligible  compared  to  the  phugoid. 
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If  a  camera  record  at  100  frames  per  second  were  made,  the  pictures  at  each 
l/lOOth  of  a  second  vould  give  eunple  data  firam.  vhich  to  plot  an  erperimentcQ. 
curve  similar  to  Figure  |1.  Also,  since  the  slower  decaying  of  the  two  short 
I>eriod  roots  was  reduced  to  1^  of  its  initial  value  at  .7  seconds,  only  the 
first  7/lOtha  of  a  second  would  be  used  in  the  graphical  analysis  (Reference 
11)  to  obtain  the  frequency  response  of  frequencies  associated  with  the  short 
period  Bu>de. 
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COHCLDSIONS  . 

The  uae   of  harnessed  models  in  circling  flight  for  evaluating 
longitudinal  stability  of  full-scale  airplanes  has  been  investigated  quali- 
tatively. 

It  was  determined  that  this  method  is  Inadvisable  for  finding  elev- 
ator power  and  stick-fixed  neutral  jKsint  because  of  the  variations  required 
in  test  power  settings  and  because  of  continuoxis  slipstream  disturbance. 

The  method  appeared  favorable  for  obtaining  transient  pitch-angle 
response  to  a  pulse  elevator  li:Q)ut,  employing  a  high-speed  motion-picture 
camera  to  record  transients. 

There  is  a  definite  lower  limit  on  the  size  of  models  which  can  be 
used  in  this  method,  considering  dynamic  similarity,  due  to  the  critical 
weight  requirement  nn^   the  necessary  mechanical  installations. 

The  use  of  tin  indoor  flying  site  featxiring  a  raised  platform  for 
the  operator  is  reccamaended.  This  would  permit  Initiation  of  transients  at 
the  datum  plane  of  A^=0  (flying  wires  horizontal)  and  would  prevent  destruc- 
tion of  the  model  in  the  event  of  ensuing  iinstable  phugold  oscillations. 

Previous  to  recording  data  the  model  should  be  flown  below  datum. 
To  preclude  effects  of  waJce  disturbance,  data  would  have  to  be  recorded  with- 
in one  lap  of  datum-level  flight. 

If  further  investigation  of  this  technique  is  Intended  with  the 
Cessna  configuration,  it  Is  recommended  that  a  model  of  at  least  twice  the 
size,  i.e.,  one-fifth  scale,  be  used. 
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APPSHDIX  A 

Development  of  general  equations  of  motion  to  include  harnessing  forces. 
This  development  is  essentially  similar  to  that  of  Reference  ^^  and  uses  the 
same  asstooptions. 

Assumptions  are  as  follovs. 

(1)  I3ifi  lines  extend  in  a  straight  line  from  the  operator  to  the 
model. 

(2)  The  height  at  which  the  model  flies  to  maintain  the  harness- 
ing lines  in  the  horizontal  plane  is  referred  to  as  datxmi  height  and 
the  model  flying  at  this  height  is  flying  in  the  datum  plane.  In 
this  plane  the  vings  of  the  model  are  assumed  level. 

(3)  The  bank  angle  of  the  model  (a^)  because  of  the  action  of 
the  harnessing  wire  guide  on  the  inboard  ving  tip  equals  the  azigle 
the  lines  make  with  the  horizontal. 


(k)     For  small  displacements  of  the  model  from  datum,  s"»n^W  -^*, 
cosii^  r  1,  and  F  remains  constant. 

(5)  The  equations  of  motion  use  stability  (vind)  axes. 

(6)  Aerodynamic  effects  due  to  cxnrvatiire  of  the  flight  i)ath 
are  neglected. 
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The  equations  of  longitudinal  motion  are  as  follovs: 

21F^    =  w^^/  =  -  Co'i/9V*S   -VVsirv*    ^  Ct  V^/^VS  (1-a) 

^Fe.    = -waVJt  =  -  Cu  ^pV\5   +  Wco&dT  -^    AvvW'^siwA^  (1-b) 

'  R 

S-C^^     -    l,<(ck-Hii)    =    Cm  ^/j^pV-JSc  (1-c) 

l^iese  equations  are  tlien  linear! ted  In  the  usual  fashion  by  considering 

each  tena:   F  (v,  y,  oc,  (J),  •  •  •  )  =  ^£.  A.V  -*-  ^  zi^ot  -4-  ^  a.  y  -i-  "b£  ^(^tV 

t>V  be*.  b^  2)^4 

After  letting       Cq  =  Co  -Cr  this  results  in 

(2-h) 

'C^^^y^py^S<i^^   -  e».^^5l^V''5c^.O  -C».^^V^pV!Sc^k>^  =  0   (2.c) 
In  order  to  put  the  term  A^  in  terms  of  the  longitudinal  variables, 
substitute  ^^  -  s  (vhere  V^  is  the  vertical  distance  from  datinx) ,  and  dif- 
ferentiate vith  re8X>ect  to  time.  Neglecting  higher-order  terms  this  gives: 


»  R  dL'X 

However,  since  initial  conditions  specify  level  flight  at  dattaa  level, 
TJ'o  -  O  ,  and    /.  4l^    -    Ay     ,  and  from  equation  (l-a),  -Vl6\v\^<^  '%S^> 
therefore  eqxoatlon  (2-b)jL  becomes 

WnVv    ^    C/<^^  Ay  -  CupVSAV  -Cc^<^Sa^-»'    ^"^^    "° 
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Dividing  through  equation  a  and  b  by    /^pV^^  r^S»    and  equation  c  by 
^  pV*"5c  =  %5c.  and  making  substitutions  as  follows: 

and  using  the  operator  ^  to  indicate  differentiation  vith  respect  to  ^'V* 
and  considering  all  «,'*'     as  being  the  displacements  from  initial  conditions 
(i.e.^  ^<*,  ^^    )  gives  the  linearized  non-dimensional  form  of  the  equation: 
niese  are  e^ressed  in  operator  form. 


V     -    o 


(3-a) 
(3-b) 


In  the  analysis  following,  Co^  and  C  »*^^  will  be  considered  small  and 
neglected. 
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AFPSMDIX  B  #- 

Derelopment  of  characteristic  equation  and  tranefer  function  ^ 
Solution  to  the  equations  of  zootion  (Equations  3^  Appendix  A)  is  nade 
by  assijming  a  solution  of  form 

vith  the  ■^^^•i^*  =  o   for  the  stick-fixed  condition  with  i>«  measured  from 
equilibrium  elevator  position. 

^Ihis  results  in 

(Co*  ^  \)>x.    -v  iCCdoc  -CcW.  ^  -L  Cu  9.  =0  (^-a) 

C^Axx.  +  [X%(U%o^^X^-P]c<,-[A\cSx^^]e, --0  (i^-b) 

For  these  honiogeneous  algebraic  equations  to  have  other  than  a  trivial 
solution  the  determinants  of  their  coefficients  must  vanish. 

Expcmaion  of  this  Srd-order  detenoinant  gives  an  equation  of  the  form 

b^\^  V  \>X   V  h^\'  ^   b.X^  ^  bA  ^V5,  X  O  (5) 

By  defining  the  combination   ^^A»  •»•  Cv^^j^.^  and  Cu^  +  Co"-  C«^  , 
the  coefficients  of  the  characteristic  equation  above  are  expressed  as 

b^  -  -  I 

b^  --  n  -  5.C0*  -  Cu^ 


^       «-        H       -"       »»  •-         Vx 
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Solution  for  the  traiMfer  function  S.  to  a  step  input  Is  made  by 
using  the  Laplace  transfona  on  the  equations  of  motion  (Equations  k). 
Since  all  initial  conditions  are  assumed  zero  for  the  variables  used,  and 
the  transform  of  a  step  function  is  unity  in  the  Laplace  tx^uisform  as  de- 
fined in  Reference  3»  this  results  in  equations  6  below. 

{CD->r>)M.    V  l-(Coo,-C.W  +   i-Cu©  -o  (6-a) 

2.  1. 

CuXju  -^[A^-wfco^o^yN^-rjoc  -[A%CDA^-f]e  =o        (6-b) 

Bote:  Although  the  left  sides  of  the  eqxoations  above  appear  similar  to 
equations  k,   in  equations  6,   u.,o(.  and  9  are  variables  of  /yt   vhile  in 
equations  k,   -^i^^.,  andd|  are  constants.  Solving  equations  6  by  determin- 
ants for  ^  gives  a  solution  of  the  form: 

^{)^^   Q-jA^  ^   <Xv  A""  -^^  <x.\  •»-  g, (7) 

**«     b^X^  ^  b^A"  ^  baX'  V  Vd.K"  ^b,X  ^  b. 

Where  the  b  coefficients  of  the  denominator  are  those  of  the  character- 
istic eqiiation  and 
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APPEMJIX  C 

Since  the  characteristic  eqiiation  is  the  denominator  of  the  transfer 

fimction  ^(X)  /  all  solutions  vill  be  made  for  the  transfer  function  only, 
o 

All  solutions  vill  be  based  upon  the  following  so-called  standard  flight 

conditions  (all  angles  in  radians): 

V   »  58.7  ft/sec   (  ^o  n^jh) 

Cu  =    .232   from  y± 

Cp  a    .0616  from  vind  tunnel  polar  (Reference  8)  corrected  by 
addition  of  vire  drag  (Appendix  D) 

/A  •  25.6 

^  =  .081^3 

*VN  B  .0^4.76  slug  ft2 

C»^^«  -.4o6  from  Reference  8 

^  »  .869  from  Reference  8 

^0^-    e-rrf^Xi.'      .116 
The  following  various  thrust  conditions  were  considered  using  harnessed 
flight  equations: 

Condition  I  "  C-r  "  O      constant 

•Condition  II  -Ct-C©     at  <xl\oc   and  ax   (Cd^=  o| 

Condition  III  -  Cv  -  C^     at  equilibrium  o<^  emd  aa.  ,  and  is 

constant  at  that  value  ( c  d^  -  C  o^* 


W     =     1.53  lb 

R  s 

37.5  ft 

S     •     1.59  ft2 

r  B 

.214 

b     •     3.28  ft 

1  • 

.119 

AR  s     6.76 

p- 

.00238 

C     «       .1^92  ft 

^^.• 

4.62 

PRINCETON  UNIVERSITY  PAGE         36 

AERONAUTICAL  ENGINEERING  LABORATORY  report    181 


These  conditions  prodixce  the  foUovlng  changes  in  the  cosrponents  of 
o^i,    and    h-^  : 


Condition  I 

Condition  II 

Condition  III 

^* 

1.0 

1.28 

1.28 

Co** 

.0616 

0 

0 

Co: 

.116 

0 

.116 

2. 

.0277 

.0269 

.0269 

Cwvj 

-.751 

Cv*a« 

-.1^9 

C*^tlo^ 

-.Olt65 

The  folloving  tenoB  depend  upon  tail  efficiency   y\^  .     Ihe  Condition  values 
vere  obtained  frcn  Beferencee  7  s^&d  8  and  corrected  for   i\f 

-.190  -.190 

-.059'*  -.059** 

Solution  of  the  transfer  functions  for  harnessed  flight  gave 

Condition   I:  ^(M  -  -Q»91  (A^-*- 2.37  ^  ^_     .2^2  X  -^  .00^) 

4>   '   (A'  V  6.99  N**  ^  8.2if9\*  ^  .7559  V-  V  .5655  ^  *  .018) 

Condition  II:  t  (^A  ^  -ll.'»3  (A""-*-  2.31  ^  ■*;  .i;^37) 

A     (  A**  •*■  7.577  A^*  +10.165  K-  -^r   .U125  /\  -».  .702) 

Condition  1X1:  ^(>C\  =  -11.^3  (A^-*-  2.31^  '«'  .l'»17) 

&     (V  ^  7.577A*->-  10.152  A^  ^   .3823  A  -^  .702) 

Factored  forms  of  the  preceding  are: 

Condition   I:  ^(X\-.    -8. 91  (A'*-  .0262)(Av  .08i»)(\->^  2.2^6) 

*  ^'   (k-.0327)(Avl.l3S)(\Z5:53)(Vi'.002\-t.2627^) 

Condition  II:  & /x\  ^  -ll.i»3  (  X  ^   .065)(A  ■>-  2.24^) 

6^  '   (X^  1.705)(A+ 5.884)(K^-.012  K  ■*-  .26i^) 

Condition  III:  tW=  Same  as  Condition  II 
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Condition  II  vas  established  as  the  most  proper  condition  and  the  trans- 
fer function  vae  re-solved  for  iinhamessed  flight  by  cofiaputing  the  coefficients 
f  or  ^  -  O  . 

Condition  11*-     0  ( x\   -      -11.1^3  (A-*-   .012)  (A  v  2,298)        

Free  flight:         6       '  (A  v  1.70l*)(/\  +   5.875)(A"  -  .OOISA  ^    .115^} 
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APPBHDIX  D 

Calctdation  of  drag  effects  iqpon  the  harnessing  wires.  This  solution 
assumes  no  outward  yaw  of  the  model.  Axes  are  as  shown. 


Ay 


H 


-»rt>-: 


cab 


Z> 


T  -  line  tension  •  F  = 


Cj)  «  1.2  two-dimensional  flow  Cj)  of  cylinder  at  low  Hg  based 
on  frontal  area. 

Sw  •  area  per  foot  of  wire  normal  to  velocity 

'     .0007  for  .008"  wire 

=  .001  for  .012"  wire 

^   •  l/2p'^r''  where  at  is  velocity  of  wire  at  station  x 

AT*  7 
B 

Coxislder  an  eLament  of  wire,  length  cAx  ,  loaded  "by  the  drag  Kx'-cLk  . 


Tcos 


(« ^  t) 


&\V\ 


^•^^a!) 
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Taking  the  summation  of  forces  in  the  drag  direction,  aesimiing  T  reanaixB 
constant  over  the  element. 


c\x  I 

T  i\v\  O     -Kx^-Xsivve<-os-  d&    -  Tcos  d  s  i"v\  4_®     -  o 

dlx.  dix 

cix  dx.  dx  cAx*- 


therefore       ^     ^ 

r^  X      -   T  e\  1 


.^      -  O 


Integrating  gives 


^  ''"    V  C  x  ^  D 


dv^     ^     -     V£_x^     ^    C 

ax  "i-r 

Integrating  again  gives 

^_ 

Solving  by  constants  for  end  conditions: 

C  --    'lA' 
Thearefore  i^"^ 

^1.  =  .     --331 

?  Vw 

For  tension  shared  by  all  vires  equally  (most  conservative  condition) 
©^^^  =   (-13. 72)(. 00234)  57.3  •  1.8»fO 
Drag  of  vires: 

Operator  v^  JL     '  a 


Ik  -"-o 


O <=>- c 


Control  handle 
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"iTS 


I^>M,  -   Ji2i-001)j^K^l  «  .0619  lb  for  .012"  vlre 

-  (3)  ^^->.  »  *lhk  lb. 

Determination  of  proportion  of  wire  drag  adding  to  airplane  drag: 
Centroid  of  vire  drag: 


Xo   = 


Therefore         .^  m  .  n.r 

3 


^  2x 


H    Us 


Therefore  the  vire  drag  centroid  is  located  3A  °^  "^^  distance  along 
the  radius  from  the  operator. 

D  vires  •  3A  (.1^)  ■  -loa^ 

Converting  to  conventional  coefficient: 

Co^  «  008  =  .0166 


Total  model  drag  coefficient: 


s  .6ko  ^   .005  ^  .0166  s  .0616 


To't^l  <iva.tt: 


D  =  Ct>  «:i^S  -    .Ottt(t.S2.^   -  .HOI. 
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Calculation  of  Line  Tension: 

Centrlfxogal  Line  Force     T,     =    ^^^     -     ^  )L 

•  .0476     .  (58.7)2  •  i^.38  lbs 
37.5 

Side  Force     ®  S*   o^  fl  -^^  =  ^^(»P     *  'OOS  x  5     =  .O^fO 

Y,   -  C^c^jS    z  (.01^0)6.52  =  .26  lbs 

Total  line  Tension  ^  k.6k  lbs. 


PRINCETON  UNIVERSITY  PAGE         42 

AERONAUTICAL  ENGINEERING  LABORATORY  report    I81 


APPEBDIX  E 


Frequency  response  calculation  of  model  In  harnessed  flighty  C^  -  ^o! 


Cooponente  of  the  Lw\-A»v^  plot  au*e  given  Figure  S  and  are  added  to 
obtain  the  ccanplete  frequency  response  in  Figure  10.  This  result  is  re- 
plotted  on  cartesian  coordinates  in  Figvtres  7  and  8  . 


"^  K  ^  Z.^il 


.0JZ<T^.oT^.31  -  l.Sdl^ 


Uwv  t(^u>^  ^  2.ofl<r«j.oi'i"i  *  2.oi^.,^7T/IEi  +  ^o^'^'-JTTTTE'Y 


-^°^-RS  -  -H.r;(g^-..o.../[.(^jj.(„,.v 


(()  e^  —  C^**^^  "  X^ol^~*  <^   ^  JLa^-'  CO    _  i:^)L^~'_co_ 
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APPENDIX  F 


TranBlent  response  by  Heavlslde  method: 
Harnessed  flight,  C,  =  Co  : 

e  _    -xiM^^ ^  2.31  ^  ^  .\w^9\       

6    '     [^x"   •*-    7-577  A*  ^  10.165 /.^  ^   A125A  -^  .702] 


P(o\  -    -ll.U3(.ll»59)     -     -2.376 
Q(o>  '  .702 


A.  =    -1.075: 

P(A.^    =     2.907  -  3.939    +   .1^59 
=  -.886 

Q'(\V  l^A,^^-   3(7.577)A"  ^    2(10.l65)/\,  ■»-    .iH25 
=  -19.82    -*-    66.079      -      3^.66         ♦         .1H25 

'<^^^-^     =         .e86(ll.i^3)  «     -.78M+ 
/\,Q'(N,)  -1.075(12.01) 

^x  -    -5.885: 

PCAv^  =    34.62    -    13.592     ■»■    .1^59 

=  21.17^ 

^CA,'^^   -814.8+786.95  -    119.96   ^    .i^l25 
-  -147.40 
^^0^^  _    -11.43(21.17^)  r     .2790 


^--i^'  -    -11.43(21.17^) 
A,.0>*(A.V  -5.884(147.40) 


|(^^  -     -2.376  -    .7844e^-°'^5*/'r^.2790e5-88^*^V 
o 

(Short  period  mode  only) 


I 
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Long  period  mode  of  harnessed  flight,    C-^  =  C  o      : 

f(*:Vi^.l2fe*°^'^^^ivv(7it  ^    135.5°) 

o 

The  coanplete  solution  for  this  condition  is  plotted  in  Figure  1 1 . 

Short  period  transient,  free  flight,  Cr  =^5: 

|(t^--  -2.383  -  .559e"^-'^0^*''^-  .2826e"5-Q'^5*V 
o 
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APPSRDIX  G 
Sample  calciilatlon  of  Moment  of  Ii^rtia   l>f 

NN  --  Model  weight  "-  1.53  lb. 

I  -  Pendulum  length  -  1.23  ft. 

^  =  Swinging  gecir  weight  -    .0331  lb. 

Si    -  Swinging  gear  i>endulum  length  -  1.053  ft. 

T  =  Measured  period  of  oscillation  =  1.335  eec. 

V  -  Model  volume  -  .216  ft.^ 

p'  Air  density  =  .00238  slugs/ft. ^ 

5=  Gravitational  acceleration  =  32.2  ft 

sec. 2 

I^  -  Inertia  of  swinging  gear 

-  riJL±.     -       .000975 

where  t.  -  measured  period  of  oscillation  of  the  gear  alone  ■=■   l.jX  sec. 

Substituting  in  the  above  eq\iation  for  l^   : 

I>j  =  ,012k   slug  ft.^ 


Fv* 


O^M.V-g      C 


"lO-llf* 


♦  f 


M     W    tWTIC 


n 


AHUi- 

•M«  (mc      fo$Fi.»^<  ) 

IS*  j«    so  n 

ST.tn.iieN 

i«  «<••     - 

riN 

t  «M        - 

El.CV*TOt     (»LUS   TMI 

«    •• 

tuioct 

»  74?       - 

«iLetoa  1  TOTAL) 

14    0* 

rtAft 

t.'M     •• 

CLHATM    TU 

«*r     - 

ANGLES  — 

mClOEMCC     tTAtlLIIl*    -t^° 

iMCiocaci    «riNu 

♦  1° 

DMWOIIAk 

«  1  ' 

tueott  TtAvcL 

t  It  ••• 

u.twKTo«  nimL 

I    JQ     jl 

AiURoii  riUvci 

♦   «f  •  •  1  • 

-  i4  ■  il  ■ 

rukP  iptn 

-  «0'  -•  f  ■ 

(UvrnBTTAt 

IS-  Il 

MMk            eti 

tfJl                  ti 

•JA            tirs 

Mti.tM.          rt"(M*ii 

■EUUlU 

f«PTT          ir« 

taoM         i«so 

I 


8l<* 


»-\ 


-«    a   >- 


(,   *(*«    SMt 


T 


t  n»o«T  »p»ii 


-»5ll"—  -     -■-*' 


tsi 


I' 


M'-tC' 


. 

H- 

\     i 

1 — 

1 

»<tL« 


'■4000< 


>-. 


M     1  I      -    I  >     I 


If' 


I 


Figs.  1  i  2 


FIGURE  1.  -  DISASSEMBLED  MODEL 


FIGURE  2.  -  COMPLETE  MODEL 


■•v^•r*J^1 


Figs.  4a  &4b 


BRAKE  ("OFF"  POSITION) 


F'iEL  REOULATOR 


PUSH  ROD  RETURJ!  SPRING 


FIGUHES  4a,  Ub.-   TOP  VIEWS  BEFORE  ASSEMBLY  SHOWINO  INTERNAL 
MECHANISMS 


Figs.  4c  8^  4d 


FIGURE  /+c.  -  VIEW  FROM  RIGHT  SIDE  BEFORE  ASSEMBLY 


FIGURE  W.  -  SHOV/ING  PRIMARY  CONTROL  WIRES  AND  CONTROL  HAMDLE 


I 


Figs.  4e  8i.4f 


FTO'IRE  4'^.    -  SHOV/TNG  METHOD   OF  W^T)TWO  RI'BBPIR  MOTOR   OF  ELEVATOR 

PULSING  MSCHAMTSM 


L  ;UitE  ur.    -  ...ODEL  CCriHONEKTS ,    DISASSEIJ3L='n  v;  '  .. 
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DETAILS    OF    Hi:/\5UR1WG 
APPARATUS 


hicr^'-      Tr)^^ey    wire    oiv%d 


M  5  -  voiV  bu.ttt'ru 
wi  tw    cent e.  r   ttvp 

I 


If        f 


r— ^ 


IT  Z^.bM 

FlaiK  bu.lb 


^ 


Figure  5 
Circuit    D  I  u.  >  r-  0.  r/^ 


iljh 


22.5V 


AY'<''pl'"^'er 


resistor 


IMHO-n. 


^ 


r 


Fig  are  b 
General     Arrange  orient 


^"1  y '^na  w  ires 


OpGV&Tor 


0:»Ct  i\oov^o>(^»^» 


i 


I 


Fitiuve  7 


\ 


i 


FUurc   8 


I 


o 


F»au.r€.  *1 


Loc^   Ma^r^ituiHe   -    Dec'ibe^s 


o 

I 


Figure   lO 


ni 


< 


Fiaure    It 


O^^     ff>      00        ^»  lO  IT)  •«■ 


0li^ 


1 


CD  -0 


«— I    C7>     oo     r^       <^         LO  ^ 


o>    00     r^      <A        i''  ^ 


Fi'^ure   12 


'■\    '-l"!.i'-    I'l^    "jil 


?  e'f^'^^Le.aM  &mB5.&^t  5tK(t^tir  ^fe5.4- 


35«5«£^t)] 
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